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"q ' Abstract 

^ | We study the topological signature of euclidean isometries in pair separations his- 

tograms (PSH) and elucidate some unsettled issues regarding distance correlations 
between cosmic sources in cosmic crystallography. Reducing the noise of individual 
j_j ■ PSH's using mean pair separations histograms we show how to distinguish between 

topological and statistical spikes. We report results of simulations that evince that 
topological spikes are not enough to distinguish between manifolds with the same 
set of Clifford translations in their covering groups, and that they are not the only 
signature of topology in PSH's corresponding to euclidean small universes. We also 
show how to evince the topological signature due to non-translational isometries. 
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1 Introduction 

The method of Cosmic Crystallography (CC), devised by Lehoucq et al. 0, looks for 
distance correlations between cosmic sources using pair separations histograms (PSH), i.e. 

* E-mail: german@cbpf.br 
^E-mail: reboucas@cbpf.br 

* E-mail: teixeira@cbpf.br 



1 



plots of the number of pairs of sources versus the distance (or squared distance) between 
them. These correlations arise from the isometries of the covering group of the 3-manifold 
used to model our Universe and so they provide a signature of its global spatial topology. 
In this way CC is potentially useful to investigate the shape and size of our Universe. It has 
recently been shown by Gomero et al. how to calculate the topological signature from 
these distance correlations in a very general geometrical-topological-observational setting. 
It turns out from the major result of Ref. |2| that correlations due to Clifford translations 
manifest as spikes in PSH's, whereas other isometries manifest as small deformations of 
the expected pair separations histogram (EPSH) of the corresponding universal covering 
manifold. 

The major result of Ref. || has a striking consequence for universe models with hyper- 
bolic spatial sections. Indeed, since no hyperbolic isometry is a Clifford translation, there 
are no topological spikes in PSH's corresponding to low density universe models. Thus, 
at first sight, these histograms seem to give no reliable information of the topology of 
the spatial sections in these models of the universe. The absence of spikes in PSH's from 
hyperbolic universes is by now well understood and has been confirmed by simulations 
performed by Lehoucq et al. || and Fagundes and Gausmann [Q. It remains, however, 
to understand the topological signature of hyperbolic isometries in CC. 

The implications of the results of Ref.|| for PSH's from flat universe models seem to 
be less well understood. It has been stated in |4| and || that every euclidean isometry 
which produces T-pairs in a given catalog will give rise to a spike in the corresponding 
PSH. This statement, however, is in clear contradiction with the fact that only translations 
produce spikes 0. Moreover, in studying the applicability of CC to closed flat models 
of our Universe, Fagundes and Gausmann reported a PSH for a manifold of class 
C?6,[] therein called model E4, which exhibits a significant peak at (d/L) 2 = 5. That 
paper suggests that this spike is generated by an isometry of the covering group of the 
manifold considered, and this interpretation was again suggested in Ref. H. Nevertheless, 
according to Ref.||, since there is no translation that would produce the peak at (d/L) 2 = 
5, one immediately concludes that it must be due to statistical fluctuations, and so 
it is not of topological origin. A definitive elucidation of these unsettled issues would 
be useful because it would clarify the actual signature of euclidean non-translational 
isometries in PSH's. Indeed, by performing simulations we will evince in this letter that, 
contrary to what is suggested in [| and 0, topological spikes are not the only signature 
of topology in PSH's corresponding to euclidean small universes. Besides we also show 

1 In this letter we use the notation of Rcf.p to denote families of flat compact orientable 3-manifolds. 
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through simulations that non-translational isometries do not manifest as less sharp peaks 
as suggested by Fagundes and Gausmann ||, but as broad and tiny deformations of the 
PSH corresponding to the simply connected case. Our results here are supported by and 
in agreement with the general theoretical developments of Ref . P] . 

Actually, the major purpose of this letter is to show how to use the MPSH technique 
described in [|] for studying the topological signature of isometries in PSH's. After a 
brief review of the techniques developed in || we first compute MPSH's for a manifold 
M of class Qq and reduce the statistical noise to a level that allows the identification of 
topological spikes. In this way (i) statistical spikes that may be confused with topological 
spikes are removed, and (ii) topological spikes that are masked by statistical fluctuations 
in individual PSH's show up even when there are few T-pairs corresponding to them. 
Incidentally, point (i) makes clear that actually there is no topological spike at (d/L) 2 = 5. 
As an additional application we construct an EPSH for the minimal 3-torus that covers 
M and plot the difference between this EPSH and an MPSH of M. Since the covering 
groups of this 3-torus and that of M have the same translations, then PSH's for these two 
manifolds would exhibit identical spike spectra 0. So this difference yields the topological 
signature of non-translational isometries of the covering group of M plus some statistical 
noise. For comparison, we also plot the difference between an MPSH and an EPSH, 
both for the 3-torus, obtaining as a result essentially statistical noise. This indicates 
that, within the accuracy of the simulations, topological spikes are the only topological 
signature in PSH's for a 3-torus. 

2 Pair separations histograms 

Here we briefly review some results obtained in Ref. 0, and extend them to the level 
needed for the development of this work. We begin by describing what a pair separations 
histogram (PSH) is, and then show how to construct mean pair separations histograms 
(MPSH) with simulated catalogs. We end with a brief explanation of the expected pair 
separations histogram (EPSH) and its use in determining the topological signature for 
non-translational isometries. 

To build a PSH we simply evaluate a suitable one-to-one function / of the separation r 
of every pair of cosmic sources from a given catalog C, and then count the number of pairs 
for which these values f(r) lie within certain sub intervals. These subintervals are all of 
equal length and must form a partition of the interval (0, f(D)], where D is the diameter 
of the observed universe corresponding to the catalog. A PSH is just a normalized plot 
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of this counting. We will take the function / to be the square function as is usual when 
dealing with flat models in CC, in order to compare our plots with those found in the 
literature. 

It is convenient to have a formal description of the above procedure. In considering 
discrete astrophysical sources, the observable universe can be viewed as that part of 
the universal covering manifold M of the space-like section M of spacetime, causally 
connected to an image of our position since the moment of matter-radiation decoupling; 
while, given a catalog of cosmic sources, the observed universe U is that part of the 
observable universe which contains all the sources listed in the catalog. So, for instance, 
the observed universe corresponding to a catalog covering the entire sky is a ball of radius 
given by the redshift cutoff of the catalog, while for a pencil beam catalog it is a thin 
cone with vertex at an image of our position. Interesting observed universes that may be 
explored in the context of CC are thin spherical shells; these observed universes correspond 
to catalogs with approximately equal upper and lower redshift cutoffs 0. 

All the sources contained in U can be observed in principle, but due to observational 
limitations a catalog consists only of part of them. Our observational limitations can 
be formulated as selection rules which describe how the catalog arises from the set of 
observable images. These selection rules, together with the distribution law which the 
objects in M obey, will be referred to as construction rules for the catalog C. It should be 
noted that the above definition for a catalog fits in with the two basic types of catalogs 
one usually finds in practice, namely real catalogs (which arise from observations) and 
simulated catalogs, which are generated under well-defined assumptions that are posed 
to mimic some observational limitations and (or) to account for simplifying hypotheses. 

To construct a PSH one begins by dividing the interval (0, D 2 } in m equal subintervals 
of length As = D 2 /m. Each subinterval has the form 



Given a catalog C of cosmic sources and denoting by rj(s) the number of pairs of sources 
in C with squared separation s, a PSH is then obtained plotting the function 



where N is the number of sources in C. Note that with the same catalog C we may 
obtain different PSH's simply by taking different values for m. The sum in (|2.1f) is just 




i = 




and is centered at 





(2.1) 
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a counting and the coefficient of the sum is a normalization constant, so 

m 

£$( Si )As = l. (2.2) 



i=l 



Although it is usual in CC to refer to the plot of the function $(s$) as a PSH, for theoretical 
purposes it is more useful to define a PSH as a function given by ( |2.1|) . From now on we 
will always refer to $(sj) simply PSH. 

Any single PSH is plagued with statistical noise that may mask the topological signa- 
ture. The simplest and most obvious way to reduce this noise is to use the MPSH which 
is described as follows. Consider K comparable catalogs Ck (k = 1,2, ... ,K), with ap- 
proximately the same number of cosmic sources and corresponding to the same manifold 
M. Let their PSH's, for a fixed value of m, be given by 

where Nf. is the number of sources in Ck and rjk(s) is the number of pairs of sources in Ck 
with squared separation s; then, the MPSH defined by 

<*(*)> = ^E**(*i) (2-4) 
fc=i 

contains much less noise than any single PSH, and clearly contains the same topological 
information. Indeed, elementary statistics tells us that the statistical fluctuations in the 
MPSH are reduced by a factor proportional to 1/ Vk, which makes at first sight the MPSH 
very attractive. In Sec.|3| we apply this technique to discriminate between topological and 
statistical spikes in PSH's corresponding to an euclidean compact manifold. 

As shown in Ref. |§ (see also Refs. 0, 0), in the limit K — > oo the MPSH approximates 
very well to the EPSH which is an "ideal" PSH, i.e. a PSH with the statistical noise 
completely removed. Equation (4.15) of Ref. [Q] [or equivalently eq. (2.11) rederived in 
Ref. P|, wherein N = n(n — l)/2 denotes the total number of pairs of cosmic images] can 
be rewritten in the form 

<$> exp (*) = + j^-j [ (*) - (*) ] + JTTJ E v* [ &Lp ~ ®ttp ] > 

(2.5) 

where T is the covering group T of M without the identity element, N is the mean value 
of the N k , $ex P ( s «) i s the EPSH of the corresponding simply connected case, N u is the 
expected number of uncorrected pairs and v u = 2N U /N. In (|2.5|) &e Xp (si) = F u (si)/As, 
where F u (si) is the probability of an uncorrelated pair to be separated by a squared 
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distance that lies in Jj. For each covering isometry it is also defined a number u g = N g /N, 
where N g is the expected number of g-pairs in a catalog with N sources, and a distribution 
function & 9 exp {si) = F g (si)/As, where F g (si) is the probability of an observed g-pair to be 
separated by a squared distance that lies in Jj. 
Now within the approximation^ 

KIM) « . ( 2 - 6 ) 

the EPSH reads 

d> exp ( Si ) « + j^-j £ ^ [ " *Sp(*) ] • (2-7) 

The general underlying setting for performing the calculations involved in (|2.7|) is the 
assumed existence of an ensemble of catalogs comparable to a given catalog C (real or 
simulated), with the same number of sources and corresponding to the same manifold 
M. The construction rules permit the computation of probabilities and expected values 
involved in ( |2.7|) . 

Let T t C T be the subset of all Clifford translations of V (i.e. all the isometries g G T 
such that for all p G M, the distance \g(p)\ = d(p,gp) is independent of p). When g G r t 
we have 

KvM) = % » (2-8) 

where 8^ is the Kronecker delta, and i g is the position of the spike due to the translation 
g G T t , i.e. \g\ 2 G J{ g . Then one can write (|2.7f) as 

$ exp ( Si ) « + ^Lp(^) + , (2-9) 

where 

= ^7 E (2-10) 

aer t 

is the contribution of Clifford translations to the topological signature of the EPSH, and 

<p(*)= AF3T S (2-H) 
JV 1 9 er\r t 

is the topological signature associated to the non-translational isometries of T. 

It is clear from (|2.9| ) that manifolds with the same translations in their covering groups 
will exhibit the same spike spectra given by (|2.10|) , the only difference between their 
2 This approximation is justified a posteriori by Fig. 5a for the case of a 3-torus. 
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EPSH's being the topological signature associated to non-translational isometries. More- 
over, in the euclidean case, from (|2.9| ) one can always write 



with 



« W*) ~ C( s «) > ( 2 - 12 ) 



*%r(8i) « + f4p(«<) ( 2 - 13 ) 



being the EPSH of the minimal 3-torus that covers M. 

Now, using the fact that <$(sj) > ~ $ e z P (si) one gets another approximate expression 
for the topological signature of non-translational isometries of small flat universes, namely 

" • (2-14) 

This expression can easily be numerically evaluated since the MPSH can be obtained with 
computer simulations, and $^ ts (s*) is given explicitly by 



stores/ 



1 ^ 1 



V <?er t / v ' g er t 



In Sec.f| we use ( J2.14p and ( J2.15|) to evince the shape of the topological signature of 



non-translational isometries for an euclidean closed manifold. 

From now on we will consider only trivial construction rules, i.e. we will assume 
that cosmic sources are uniformly distributed in space, and all cosmic sources present in 
universe models, up to a given redshift, are recorded in catalogs. Although unrealistic, this 
assumption makes easy to illustrate the general results developed in Ref.|2| and permits 
a comparison with current literature in CC. Besides, in this case we can readily compute 
^ex P ( s i) an d the coefficients v g . 



Indeed, from Ref.|T0[ (see also Ref. |Tl[) one can easily calculate for the case 

of trivial construction rules and a ball of radius R as an observed universe. For flat models 
one obtains 



Kl P (s) = |^ (2 - (4 + ^) 6(2 - , (2.16) 

where is the Heaviside function. Correspondingly, the coefficients v g can be calculated 
by simple geometrical arguments. Indeed, let U be the observed universe, i.e. a ball 
of radius R centered at our position. The isometry g transforms isometrically the ball 
g~ l (U) into the ball U, so only the sources in g _1 (W) DU have a ^-partner in U, and form 
g-pairs. Thus we have 

= VoKg-\U)nU) 

Vol(W) • 1 j 
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A simple calculation yields 

— iftB© 3 ' (2 - i8) 

where d g is the distance from the center of the observed universe to its image by the 
isometry g. 



3 Topological and statistical spikes 

Let us now show how to use the MPSH technique to discriminate between topological 
and statistical spikes by working out two examples of models of the universe reported in 
Ref. [[[J. In order to make a comparison with the plots of the upper part of fig.l in Ref.0, 
we took a manifold M of type Qq with covering group T generated by 

<*(x,y,z) = (x + l,-y,-z), 

P(x,y,z) = (-x,z + l,y), (3.1) 
S(x,y,z) = (-x,z,y + l), 

where (x,y,z) e R 3 . A fundamental polyhedron for M and a detailed construction of 
this manifold is given in Ref. [|T2f (see also [|I3| ). 



We have performed simulations for two observed universes with radii R\ = y/2 and 
i?2 = 5/6, respectively As in Ref.j7|, for each simulation in the first case (Ri = y/2) we 
put 20 objects uniformly distributed in the FP, while in the case i?2 = 5/6 we put 101 
objects inside it. In both cases we end up with catalogs of approximately 240 sources. A 
PSH for one catalog for each case is shown in fig.l, where we have subdivided the intervals 
(0,4R?] in bins of width 0.01 as in Ref. [^]. A direct comparison with the plots in Ref. 
can be done by considering unit lengths of 4200/i _1 Mpc for the first case and 7200h~ 1 Mpc 
for the second. An important difference between our plots and those in Ref. |7J is that 
ours are normalized histograms as described in Sec.§. We can see that the PSH in fig. la 
presents no apparent spike at s — 5, suggesting that the corresponding spike in Ref. is 
of statistical origin. In order to elucidate this and related issues we shall now discuss the 
MPSH's corresponding to these PSH's. 

Figure 2 shows one MPSH built with 50 computer generated comparable catalogs 
with approximately 240 sources for each case of study. It can be seen that, by taking 
the mean over 50 catalogs, the statistical noise has been considerably reduced. Moreover 
it became apparent from such plots that there is no topological spike at s = 5 for the 
first model of the universe, whereas for the second there is a small spike at s = 2 that is 
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masked by statistical fluctuations in the individual PSH's of fig. lb of this work and fig. 1 
of Ref. 0. Both MPSH's show the agreement between our simulations and the theoretical 
results of 0. Besides, it can be shown that there are no topological spikes at any odd 
integer position in PSH's corresponding to models of the universe which have M as spatial 
sections of spacetime (see Appendix A for details). 



4 Topological signature of non-translational isome- 
tries 

Let us now consider the group T t generated by the translations a = a 2 , b = (3 2 and 
c = . This group is formed by all the translations in V and is the covering group 
of the minimal 3-torus that covers M. A fundamental polyhedron for this 3-torus is a 
parallelepiped of height 2 and square base of side a/2. As it has been shown in Sec.|], PSH's 
corresponding to models of the universe that have the same translations in their covering 
groups exhibit the same spike spectra of topological origin. This result is illustrated in 
fig. 3 where it is shown a PSH, and an MPSH computed with 50 simulated catalogs in 
this 3-torus . All catalogs used here are comparable to those used in the simulations 
performed in Sec.|3| for an observed universe of radius Ri, and have approximately 240 
sources. We can note that there is no relevant difference between plots corresponding to 
M and its minimal covering 3-torus even after reducing the noise by calculating the mean 
over 50 catalogs. Thus, looking at any of these plots one cannot say whether it comes 
from a catalog of sources in a Qq universe or in its 3-torus minimal covering: topological 
spikes are not enough to distinguish between these two flat manifolds because they have 
the same set of Clifford translations in their covering groups. 

To evince the topological signature of non-translational isometries of M in CC, and 
so be able to distinguish between M and any other manifold, we make use of fl2.14| ) 



and ( |2.15|) . Fig.4 shows the graph of $g°™ s (s.;) given by eq. (|2.15|) together with ( |2.16|) 



for the case of an observed universe of radius R± = \/2. The comparison between the 
EPSH of fig.4 and the MPSH shown in fig.3b makes apparent the strength of the MPSH 
procedure. In both fig.4 and fig. 5 we have reduced the number of sources from 240 
to 120 to enhance the amplitude of the topological signature. In these figures we have 
also incremented the width of the bins from 0.01 to 0.02 just to keep the amplitude of 
the spikes as in the previous simulations. Figs. 5a and 5b are plots of eq. ( |2.14| ). The 
first plot corresponds to eq.( |2.14j ) with an MPSH, < $(s«) > , for the 3-torus while the 
second corresponds to eq.( |2.14] ) with an MPSH for M. Both MPSH's were built with 
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5000 catalogs. Since the covering groups of this 3-torus and that of M have the same 
translations, PSH's for these two manifolds would exhibit identical spike spectra; so the 
differences given by eq. (|2 . 1 4j j yield the topological signature of eventual non-translational 
isometries plus some statistical noise. Fig. 5a is essentially noise, as expected, since there 
are no non-translational isometries in the covering group of the 3-torus; while fig. 5b 
evinces the topological signature of non-translational isometries of T (the covering group 
of M). In fig. 5b, the amplitude of the topological signature of non-translational isometries 
is roughly one order of magnitude smaller than the amplitude and the statistical noise of 
the PSH's in fig.l, so this signature could hardly be eventually extracted from a single 
PSH constructed with a real catalog. We have been able to evince this tiny topological 
signature in our simulations because we have succeeded in reducing the statistical noise 
by two orders of magnitude using the MPSH technique. Incidentally, if there is any other 
component to the topological signature in PSH's of 3-tori (i.e. other than spikes), its 
amplitude must be at least three orders of magnitude smaller than that of a PSH, as 
shown in fig. 5a. 

One can obtain a better understanding of the structural features exhibited in fig. 5b, in 
particular the exact locations of the "jumps", by noting that, since V acts freely and dis- 
continuously on euclidean space, there is a minimum non-null pair separation |g| m m among 
all g-pairs, for each g e T. Thus, the terms 3>f (sj), which appear in the definition (|2.11 ) 
of the topological signature of non-translational isometries, are zero until the i-th bin 
given by the condition \g^ min G Jj. As reported in || (see fig. 2a and fig.2b therein), these 
terms start with a jump at this i-th bin, therefore one expects that the topological sig- 
nature of non-translational isometries presents such jumps at well defined positions that 
we will calculate in what follows. 

Note that all non-translational isometries of T can be grouped in the following three 
categories according to whether the parameters I and m are even (E) or odd (O) (see 
Appendix A) 

Type OE I = 2k x + 1 and m = 2k 2 , 
Type EO / = 2*i and m = 2k 2 + 1 , and 

Type OO I = 2* x + 1 and m = 2k 2 + 1 , 

where ki and k 2 are integers. 

Straightforward calculations yield 

|p(2*i + l,2* 2 ,n)|Ln = (2* 1 + 1) 2 , 
|/i(2*i + l,2* 2 ,n)|L n = (2* 1 + 1) 2 , 
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Table 1: The first five values for the squared minimum distance \g\ min of g-pairs, and 
specifications of the type of isometries together with the values of the corresponding 
parameters which give rise to them. 



for type OE isometries; 

10(2*1,2*2 + 1,71) 

|/j(2*i,2* 2 + l,n) 
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for type EO isometries; and 
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2(* 2 + ^ + l/2) 2 

2(* 2 + 77 + 1/2) 2 



2(77 -1/2) 2 , 
2(77+ 1/2) 2 . 



for type 00 isometries. Here p(l,m,n) and fi(l,m,n) are given by equation (A.l). 

Now it can easily be seen that the non-translational contributions to the topological 
signature sharply start at Si = 0.5, 1, 4.5, 9, 12.5 . . . (see table |l] for clarification). An ob- 
servable universe for a full sky covering survey and radius R\ allows non-null contributions 
from non-translational isometries with |<7|^„ < 4i? 2 = 8, and so its topological signature 
must present "jumps" at s» = 0.5, 1 and 4.5 only, in agreement with fig. 5b. 

To close this section we emphasize that the MPSH technique is a suitable approach 
to obtain the topological signatures only when one is dealing with simulated catalogs. 
Note however that the ultimate step in most of such statistical approaches to extract the 
topological signature is the comparison of the graphs (signature) obtained from simulated 
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catalogs against similar graphs generated from real catalogs. To do so one clearly has to 
have the simulated patterns of the topological signatures of the manifolds, which can be 
achieved by the MPSH approach discussed in this work. 

5 Final remarks 

In this letter we have studied, through simulations concerning closed flat models of our 
Universe, some general results obtained in Ref. 0, namely, that topological spikes in PSH's 
arise from translations, while non-translational isometries manifest as slight deformations 
of the EPSH of the corresponding simply connected case. 

We have used the technique of taking means of PSH's, described in Ref. @, to reduce 
the noise to a level that allows (i) to distinguish between topological and statistical spikes, 
and (ii) to evince the shape of the topological signature of non-translational isometries. 
In particular, we have shown explicitly that topological spikes for a cubic manifold M of 
class Qq with L = 1, appear only at even integers in plots of n(d) vs. d 2 , so the spike 
at (d/L) 2 = 5 reported in Fig.l of Ref. JF| is indeed due to statistical fluctuations, in 
agreement with the theoretical conclusion we have derived from the results of 0. 

Identifying the shape of the topological signature of non-translational isometries is 
possible only through a drastic reduction of the statistical noise followed by the elimination 
of the uncorrelated part of the EPSH together with the topological spikes. This last step 
has been attained by subtracting the EPSH of the minimal 3-torus that covers M from 
an MPSH built from M and with small enough noise. In doing so we have shown that 
the topological signature of non-translational isometries is formed by broad and tiny 
distributions, each one beginning with a "jump" at |p|^ in , for some isometry g. 

The methods employed in this work allowed the study of the nature of the topological 
signature in CC of euclidean non-translational isometries, but can be equally applied to 
any of the other two geometries of constant curvature. Note, however, that these methods 
may not be useful in applications of CC to real catalogs, since it is impossible in practice 
to construct thousands, or even hundreds, of comparable catalogs of real cosmic sources. 
For an implementation of CC along these lines for use with real catalogs, a novel technique 
for reducing statistical noise with just one catalog has to be developed. We emphasize 
however that, despite the use of MPSH's be restricted to simulated catalogs, it is a suitable 
approach for studying the topological signature of non-translational isometries in PSH's, 
and thus it is a very important tool for the understanding of the method of CC. Clearly, 
without the understanding which arises from the MPSH approach to CC it would be quite 
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difficult to establish the applicability of the crystallographic method as well as to devise 
alternative methods in a systematic way. 
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Appendix A 

In this Appendix we show explicitly that there are no topological spikes at odd integer 
positions in PSH's constructed from models using M (given by the covering group T 
generated by ( |3.1| )) as spatial sections of spacetime. From the relations [3 2 = 5 2 , Pa = 
5a = a~ l 5, and (3a' 1 = a5~ l , and from the fact that a 2 , P 2 , [35 and 5(3 are 
translations, and so they commute, one can write any covering isometry of M in one of 
the two canonical forms 

p{l,m,n) = a l (3 m {(35) n , 

(A.l) 

fi{l,m,n) = a l (3 m (5(3) n , 

where I, m and n are integers, and n > 0. For all I and m we have that p(l,m, 0) = 
p(l, m, 0), while for n^O, p{l, m, n) ^ p(l, m, n). 

The translations in V are those isometries with / and m even, regardless of the value 
of n. Explicitly we have 

p(2k 1} 2k 2 , n)(x, y, z) = (x + 2k 1} y + k 2 + 2n, z + k 2 ) , 
p{2ki, 2k 2 , n)(x, y, z) = (x + 2kx, y + k 2 , z + k 2 + 2n) , 

where k\ and k 2 are integers, and so 

\p(2k u 2k 2 , n)\ 2 = \p(2k u 2k 2 , n )\ 2 = Ak\ + k 2 2 + (k 2 + 2nf . (A.2) 



Equation ( |A.2| ) gives the position of any potential topological spike in PSH's based on 



models of the universe with M as a spatial section of spacetime. Whether a given spike 
will appear depends on the shape and size of the observed universe, i.e. on the shape and 
deepness of the region scanned by the astronomical survey used to construct the catalog. 
However, by writing Eq.( |A.2j ) in the form 



p(2k u 2k 2 , n) | 2 = 2 2k\ + {k 2 + n) 2 + n 2 , (A.3) 
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and an identical expression for |//(2A; 1 , 2k 2 , n)| 2 , one can conclude that there are no spikes 
of topological origin at odd integer positions, since \p(2k\, 2k 2 , n)\ 2 and \fi(2k 1 ,2k 2 ,n)\ 2 
take always even integer values. In particular, there is no topological spike at s — 5. This 
is in agreement with the simulations which gave rise to figs. 1 and 2, and shows that the 
spike found by Fagundes and Gausmann at (d/L) 2 = 5 is of statistical origin. 



Captions for the figures 

Figure 1. PSH's for two Qq universes. In (a) the corresponding observed universe 

is a ball of radius R\ = a/2, in (b) R 2 = 5/6. In both cases the catalogs have 
approximately 240 sources. The intervals (0, 4R 2 } have been subdivided in bins of 
width 0.01. In (a) the PSH presents no apparent spike at s = 5. 

Figure 2. MPSH's <$(sj) > built with 50 simulated comparable catalogs with approxi- 
mately 240 sources for the two Qq universes of fig.l. The statistical noise has been 
considerably reduced so that it becomes apparent that there is no topological spike 
at s = 5 in (a), whereas in (b) there is a small spike at s = 2 that is masked 
by statistical fluctuations in the PSH of fig. lb. The intervals (0,4i? 2 ] have been 
subdivided in bins of width 0.01. 

Figure 3. Part (a) is a PSH $(sj) and part (b) is a MPSH < > for the 3-torus 
universe which covers the Qq universe of figs. la and 2a. All catalogs used here are 
comparable to those used for figs. la and 2a, and have approximately 240 sources. 
The intervals (0, 4R 2 ] have been subdivided in bins of width 0.01. There is no 
relevant difference between graphs corresponding to the Qq model and its minimal 
covering 3-torus, illustrating that topological spikes are not enough to distinguish 
between these two flat manifolds. 

Figure 4. An EPSH § exp (si) given by (|2.15|) for the 3-torus whose PSH is shown 

in fig. 3. The comparison between the EPSH of the present figure with the MPSH 
of fig. 3b makes apparent the suitability and strength of the MPSH procedure. 

Figure 5. MPSH's corresponding to the topological signature of non-translational isome- 
tries given by ( |2.14f ). Part (a) corresponds to the 3-torus of fig. 3 and fig.4, while 
(b) corresponds to the Qq manifold of figs.l and 2. Both MPSH's were built with 
5000 catalogs of approximately 120 sources, and with bins of width 0.02. While (a) 
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exhibits essentially statistical noise as expected, (b) shows the topological signature 
of non-translational isometries of Qq. 



References 

[1] R. Lehoucq, M. Lachieze-Rey & J. -P. Luminet, Astron. Astrophys. 313, 339 (1996). 

[2] G.I. Gomero, A.F.F. Teixeira, M.J. Rebougas & A. Bernui, Spikes in Cosmic Crys- 
tallography, fer-qc/981 10381 (1998). 

[3] G.I. Gomero, M.J. Rebougas & A.F.F. Teixeira, A Signature for the Shape of the 
Universe, Int. J. Mod. Phys. D, in press (2000). 

[4] R. Lehoucq, J. -P. Luminet & J. -P. Uzan, Astron. Astrophys. 344, 735 (1999). 

[5] H.V. Fagundes & E. Gausmann, Cosmic Crystallography in Compact Hyperbolic 
Universes, |astro-ph/98lT368| (1998). 



[6] J. -P. Uzan, R. Lehoucq, & J.-P. Luminet, Astron. Astrophys. 351, 766 (1999). 

[7] H.V. Fagundes & E. Gausmann, Phys. Lett. A 238, 235 (1998). 

J. A. Wolf, Spaces of Constant Curvature, fifth ed., Publish or Perish Inc., Delaware 
(1984). 

[9] G.I. Gomero, M.J. Rebougas & A.F.F. Teixeira, A topological signature in cosmic 
topology, |gr-qc/991 1049| (1999). 

[10] A. Bernui & A.F.F. Teixeira, Cosmic crystallography: three multi-purpose functions, 
Mtro-ph/9904l80| (1999). 

[11] M.J. Rebougas Distinguishing Marks of Simply-Connected Universes, gr-qc/000740. 
Int. J. Mod. Phys. D 9, in press (2000). 

[12] G.I. Gomero, Fundamental Polyhedron and Glueing Data for the Sixth Euclidean 
Compact Orientable 3-manifold, preprint CBPF-NF-049/97 (1997). 

[13] A. Bernui, G.I. Gomero, M.J. Rebougas & A.F.F. Teixeira, Phys. Rev. D 57, 4699 
(1998). 



15 




Figure 1 




Figure 2 




Figure 3 




Figure 4 




Figure 5 



